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Abstract 

A collisional electromagnetic dispersion relation is derived from 
two-fluid theory for the interchange mode coupled to the Alfven, acoustic 
drift and entropy modes in a partially ionized plasma. The fundamental 
electromagnetic nature of the interchange mode is noted; coupling to the 
intermediate Alfven mode is strongly stabilizing for finite Both ion 

viscous and ion-neutral stabilization are included, and it is found that 
collisions destroy the FLR cutoff at short perpendicular wavelengths. 
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Intr oduction 

This paper is motivated by a study of the instability of the 
nighttime equatorial F region called spread F (Hudson 1974; Hudson 
and Kennel 1974a)* W^e develop a set of plasma equations to suit 
equatorial F region parameters, which are also applicable to Q 
machines and other laboratory plasmas. Although Coulomb colli- 
sions predominate in the F region, neutral collisions may determine 
the altitude threshold for spread F onset (Hudson 1974; Hudson and 
Kennel 1974a), so both have been included in our analysis. Colli- 
sional fluid equations are used (Braginskii 1965), which are appro- 
priate for parallel wavelengths exceeding the total electron mean free 
path and perpendicular wavelengths greater than the ion Larmor 
radius. Equatorial localization permits even longer parallel wave- 
lengths, so both finite electron heat conduction along the magnetic 
field (Tsai et al, 1970) and energy transfer between species must be 
included. Since this is the first of two papers on the low frequency, 
density gradient driven modes of a partially ionized coUisional plasma, 
the full dispersion relation for all such modes (interchange, entropy, 
and drift) will be derived and the modes will be decoupled here. In 
the second part of this paper we will be concerned with the structure 
of the interchange mode: its electromagnetic corrections, ion finite 

Larmor radius (FLR) stabilization (Rosenbluth, et al* 1962) and the 
relative effect of neutral and coulomb collisional damping. 

In order to treat the interchange mode properly, it is necessary 

to derive an electromagnetic dispersion relation. The interchange mode 
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; is often derived in the electrostatic approximation neglecting pertur- 
bations along the magnetic field (k„ = 0) (cf, Rosenbiuth et al. 1962). 

While this procedure gives the correct dispersion relation at = 0, 
it conceals the fundamental electromagnetic nature of the interchange 
mode. The interchange mode appears in lowest order in the ion FLR 
parameter b in the small ^ limit of the full electromagnetic 

dispersion relation; electrostatic modes such as the drift which require 
finite k^ are obtained in the large limit, where is the 

Aliv^n speed. It is incorrect to solve for the finite parallel wavelength 
(k^) corrections to the interchange mode from the electrostatic disper- 
sion relation, since the dominant stabilizing term comes from coupling 
to the electromagnetic Alfv^n mode. 

The collisionless interchange or Rayleigh- Taylor mode destabi- 
lized by gravity antiparallel to a density gradient and perpendicular to the 
magnetic field was first suggested by Dungey (1956) as a source of equa- 
torial F region irregularities. Haerendel (Balsley et al. 1972; Haerendel 
1974) included neutral collisions but neglected Coulomb collisions in the 
Rayleigh-Taylor mode, so his analysis is restricted to lower altitudes 
than typical spread F observations at finite perpendicular wavelengths 
(Hudson and Kennel 1974a). Both Dungey and Haerendel neglected FLR 
stabilization, which Rosenbiuth et al. (1962) have shown to be important 
at short perpendicular wavelengths. 

The purpose of the second part of this paper which deals exclu- 
sively with the interchange mode is, to extend the previous work in 
slab geometry to higher altitudes and answer the following questions. 

What is the effect of Coulomb collisions on the Rayleigh-Taylor growth 
rate, and how does it depend on plasma density? This is compared with 
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the growth rate dependence on neutral collisions and neutral density 
(Balsley et ah 1972; HaerendeJ 1974), What is the shortest perpen- 
dicular wavelength above the FLR cutoff, and how do collisions affect 
this FLR cutoff? This is compared with the observation that spread F 
perpendicular wavelengths can extend below the Lon Larmor radius. 

We will find that collisional particle diffusion drifts oppose the collision- 
less FLR drift of Rosenbiuth et al. (1962), and can extend the unstable 

perpendicular wavelength range down to the ion Larmor radius where 

♦ 

the approximations break down. However, the effect of both neutral 
and Coulomb collisions at long perpendicular wavelengths is to reduce 
the maximum growth rate of the interchange mode. 

Assumptions and Basic Equations 

c *>0 0 *hat a par^^ally ’oniaed slab plasma Is Inamersed in 
a uniform z-dirccted magnetic field. There is a constant vertical den- 
sity gradient in the positive x-direction and gravitational acceleration 
g in the negative x-direction. We neglect particle sources and sinks, 
zero order drifts along x and z, static electric fields and zero order 
temperature gradients. 

We will restrict our analysis to perpendicular wavelengths greater 

than the ion Larmor radius and parallel wavelengths greater than the electron 

mean free path for momentum transfer to ions and neutrals X = a /(v . + v ), 

which depends on the electron thermal speed a^ and the sum of the electron- 

ion and electron-neutral collision frequencies v . and v defined in Table 1 . 

ei en 
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Tiui following set of fluid equations then applies to the jth species, electrons 
or ions for j = e or i 
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m = m , m. = M 

e L 

Equations (1) ~ (3) are the momentum, continuity, and heat flow 
equations written in the neutral rest frame. Neutral dynamics are neglected 
for oscillation frequencies satisfying u» > V. rijn , However momentum 
and energy loss to the neutral sink are included. 

The transport coefficients for a fully ionized plasma including ion 
dynamics in an arbirtary magnetic field were computed by Braginskii (1965), 
Mis resistive, thermo- elect j*ic and electron thermal conductivity coeffi- 
cients along the magnetic field are = 0. 51, = 0. 71 and - 3, 16 

respectively, Shkarofsky (1961) has tabulated them as functions of 

for a partially ionized plasma neglecting ion dynamics, hence 
ion drift with respect to neutrals, and Schunk and Walker (1970) have 
plotted them. These exact numerical coefficients enable us to write 
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the electron equations ( I - 3 ) in the identical form to the fully ionized 
case roplacinu v . by the total electron collision frequency v = v . + v 
For example in Equation (2) we have for V. = 0 

L 

R . + R = -C m n V V 
ei en r e e e 

where C is a function of v /v 
r en ei 

Approximations in ion equations 

In the limit |k C /(O^l «1* where C = T /M is the ion acoustic 
' z s D ' s e 

speed and ~ diamagnetic drift frequency, ions primarily 

move perpendicular to B, undergoing a shear rather than a compressional 

motion (Tsai et al» 1970)« As a result, ion temperature fluctuations are 

scaled down from electron temperature fluctuations by the ion FLR factor 

b = k C /Q. <<1. We have independently checked that including ion 

i. s 1 

temperature fluctuations does not significantly affect the interchange mode 
(Hudson 1974). Therefore, we will neglect them here, thus eliminating 
the ion heat transfer equation from the set (1 - 3), 

We retain parallel ion pressure, finite ion inertia, FLR effects 
and viscosity. We use the fully ionized ion viscosity tensor (Braginskii 1965), 
since momentum exchange between ions and neutrals is treated separately 
in (4c), Our separate treatment of the two collision processes makes the 
reasonable assumption that an ion-neutral collision transfers the total ion 
mementum to the neutrals, while Coulomb collisions are a diffusion pro- 
cess in a spatial gradient. This is a good approximation for charge 
exchange collisions, and hard sphere collisions when ion and neutral 
masses are comparable, since the momentum exchange rate then 

scales as M./M 1. 

i n 
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The only contributing terms in the ion stress tensor to order 

I 

b = k^C^/n^ are (Shkarofsky _et_^. 1963) 

1 ® ^ 

V. TT= -1/(2Q.) {fp.V^ + (V P.) . V l(v X z) 

1 ’-i- 1 X X i 1-* X 

+ (zx vp.) . V V y 

J- ^ ^ _L^ 

- (ViV"D • 

- VPiV^ij • 

- l/3vj(p./v.i)y^. vj 


(5) 


The coUisionless stress term is the ion FLR effect. The first collision 

term corresponds to shear stress due to ion viscosity, and the second 

corresponds to compressional stress from the collisional relaxation of 

3 

n - n differences (Stix 1969k The compressional term is of order b 
* 1 ‘ I! 

in the ion continuity equation (Stix 1969) and will be dropped along with all 

2 

other terms higher order than b . Parallel ion viscosity like parallel ion 

heat conduction, which has been included in the collisional electrostatic 

drift-acoustic dispersion relation of Coppi and Mazzucato (1971), only 

1 

affects short parallel wavelength modes k \ '~(m/M)'^ and will be 
neglected here. The remaining ion viscosity term then determines 
the coefficient of perpendicular ion- ion collisional momentum transfer. 

We will see that this factor is reduced from unity, the coefficient of 
ion neutral collisional momentum transfer in (4c). In (5) we have 

9 

neglected coUisionless FLR terms of order b u) (Kennel and Greene 1966) 

2 

and retained those of order b v_ for application to low frequency oscil- 
lations cju/v., <<l, 

11 
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Approximations in electron equations 

We neglect electron inertia at frequencies low compared to the elec- 
tron plasn\a frequency o) < < u) . The perpendicular electron motion can be 

pe 

treated in the guiding center approximation for the modes of interest 
(k /k <<l). We therefore neglect the perpendicular pressure, 
diamagnetic drift and off diagonal heat flow terms in the electron equations 
(I - 3). It can be shown that including these does not alter the final 
result since all additional terms cancel in the perturbed equations. 

Electron thermal conductivity along the magnetic field greatly exceeds 
ion thermal conductivity in general: hence electron temperature fluctuations 
have been included along with parallel electron pressure. 

The parallel electron heat flux is 


q “ C.n Tv - C 
t e e ez x 


( n T /m 
\ e e 


^ rrT 

' "e 




Aiiain and C are functions of v / v . (Shkarofsky 1961) defined so as to 
° t X en ei 

write the electron equations (I - 3) in the fully ionized form, replacing v^. by 

V = V + V . The collisional energy transfer from electrons to ions and 
e ei en 

neutrals is given by 


= -3(m/M) (T^ - T;) = Tj (6b) 

Here we assume that ion and neutral masses and temperatures are equal, 

and that all electron energy lost to the ions is subsequently lost by ion 

neutral collisions to the neutral sink. This assumption is valid for 

V. / V.. >2 /m/ M (Hud son 1974), For v. / v.. < 2 , / m/M the ions 
in 11 in u 

prefer to give their energy back to the electrons via Coulomb collisions. 
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' becomes important when parallel wavelengths are comparable 

to the energy transfer mean free path from electrons to ions and neutrals, a 
factor of v^M/m longer tlian the momentum transfer mean free path * 

Such wavelengths are typically long compared to laboratory plasmas, so 
this term has previously been neglected in the derivation of the following 
modes (Tsai et al . 1970). However, the possibility of very long parallel 
wavelengths exists in the ionosphere, so the energy transfer term has been 
retained. 

Equilibrium 


Perkins (1973) has examined the F region equilibrium, where the 
Hall conductivity is negligible, and Coulomb collisions make no contribu- 
tion to the Pederson conductivity. In the frame where neutral winds 
vanish, the ion momentum eauation (Z) yields 
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Since v.^/Q<<! in the F region, higher order terms are dropped. 
The electron momentum equation (2) yields 

V = ^ 

B 


(7a) 


(7b) 


Perkins (1973) finds tliat only the East-West (y) component of the 
current = no ( V | j - V |c) contributes to the equilibrium. Since the 
primary electric field in the nighttime F layer is vertical due to polari- 
zation (Rishbeth 1970), the only contribution to the equilibrium comes 
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from 


A 


cT 


B 


V n X z + n e x z 
' Q 


(7c) 


which is independent of electric fields and neutral collisions. This 
is just the current due to the density gradient drift and gravitational 
guiding center drift. We will neglect electric fields in the perturbed 
equations which introduce the ^ x j[^ instability (Simon 1962) to be 
treated in a separate pap«r in the finite heat conduction limit. 

Perturbed equations 

We linearize the above set of fluid equations assuming the fluctuations 
are low frequency oscillations of the form For k L 1 the x- 

dependence of the perturbation is weak and a good estimate of the eigen- 
frequency is obtained by setting k^ = 0 (Krall 1968), 

The parallel ion momentum transfer equation yields 

7 f 1 ) g . g ■) 

zt o)ML \ u) / z cu y J 

2 

. n k n. 

+ Rrp Z S 1^ 

n (8a) 

We solve the perpendicular ion momentum transfer equation iteratively 

for low oscillation frequencies and low ion- ion and ion- neutral collision 

frequencies v.. and v, defined in Table 1 assuming v. ./Q. and 

^ ii in ® i 11 i 

b < < 1 appropriate for the F- region where ions are magnetized. 

We obtain 


y» 


V , ta C 
X I 


■ '(«/"iKi - 9/f^i "i/" 

i 0-3b^v,,[(ie^kyTjty 


(8b) 


(8c) 



where 


U) « (D + I V 


LI. 


v^/b = Vj„+0.3b VijRT. ‘^T ” h/^e 

V /b is the frequency which scales the perpendicular ion diffusion rate. 

± 

Vj.^and always appear in this set of equations in this combination. Ion- 
ion collisions only affect transport in the presence of a spatial gradient, 
where ions in a denser region encounter more collisions with other ions 
and diffuse toward less dense regions. This is in contrast to a homogeneous 
plasma, where the center of mass remains stationary in an ion- ion collision 


and there is effectively no diffusion. It further contrasts with ion neutral 
collisions where we assume that the total ion momentum is lost. Even if 
only a fraction of the ion momentum is lost in a single ion-neutral collision, 
the ions will random-walk and diffuse through the neutrals, which act as a 
momentum sink. 

For kyLr »1 the important scale length is the perpendicular wave- 
length. The different b dependences of and in v suggest that the two 
collision processes will be significant over different perpendicular wave- 
length ranges. We will subsequently find that ion-neutral collisions damp 
the drift mode at long perpendicular wavelengths. Ion viscosity will also 
be found to increase the collisional damping of the interchange mode at 
short perpendicular wavelengths. 

Substituting v^., v^, and v^. into the perturbed ion continuity 

2 

equation, assuming k L >> 1 and dropping terms of higher order than b 

y 1 

eliminates the collisional terms in (8c) and yields the following relation 
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between ion density and electric field fluctuations 

= ^ (e/STj«,o - bu, - 

+ Vy/'”)^z 


V ./n 

ei e 


(9) 

The electron momentum transfer equation yields, assuming <«/0 » 

and V /O << 1 
en e 


V = cJ£) /B 
xe y 


(10a) 


V 

ye 



(10b) , 



2 2 

where v = k a /2C v and v = v . + v . v scales inversely with the 
II z e re e e i en n 

sum of the collision frequencies. It must be small compared to the electron 
transit frequency over a parallel wavelength, k a , since k a /v = k X «1 
By assumption many collisions must occur over a parallel wavelength for this 
t'.ul 1 i uional fluid analysis to apply. 
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We substitute and into the electron heat flow equation and 

retain energy transfer from electrons to ions and neutrals. At very Low 

oscillation frequencies + ^en J collision frequencies adjust to 

the electron density and temperature fluctuations; therefore we must perturb 

the electron density and temperature dependences of the collision frequencies. 

-3/2 . 1/2 

From Table 1 v^. , while cc n^T^ is independent of electron 

density. The perturbed electron heat flow equation yields 


+ i2(l -R^j{m/M)(2v^, 

/(w+iviix) (11a) 

where 

X- 2/3[c^C^ + (l+cj2]+i2(m'M)vg^S 

+ i . (l - RT)(m/M)(vg^- 3vgi)/v|| 

(11b) 

Substituting (lOa)f (10c) and (11) into the electron continuity equation yields 
the following relation between electron density and electric field fluctuations 
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n J rj)[(u) + 1 VII )(«> ■*• iviix) * 

•" ■"‘‘y^'y/'") 

. (u) + I v\ix ) + 


(12a) 


where 


U)< 


k V 

y ey 


- kyC^nt, 


(12b) 


5 » 2/3(w02^2(i^C^)(,.R^) 


(12c) 


■ (m/M)(2vg, + Vg^)y'vii 

~X and ^ are combinations of terms which arise from finite parallel 

electron heat conduction which permits electron temperature fluctuations, 

and energy transfer between species. They are discussed further and 

plotted as functions of v / v . by Hudson and Kennel (1974b). 

on 6 1 

Electromagnetic Dispersion Relation 

To the set of five fluid equations we add Faraday’s law 


1 SB _ « 

and Ampere's law, neglecting the displacement current for ir < 

= lUj 

c ~ 


(13a) 


Vx B 


(13b) 



15 . 


* No).k;ctin^ tlio displacement current is equivalent to assuming quasi 
neutrality in place of Poisson’s equation. 

Following Mikhailovskii and Rudakov (1963) and Coroniti and 


Kennel (1970) we will neglect the parallel component of the oscilla- 
tory magnetic field. This eliminates the fast, or magnetosonie mode 
of Iiydromagnetics, which does not couple significantly to the modes 
of interest here, with perpendicular phase velocities the order of the 
diamagnetic drift. However, the intermediate and slow modes are 
well known to couple to drift and flute modes (Kennel and Greene 1966) 
so tlic AUv6n and acoustic terms will be included in our initial set of 
equati ons. 


~ ^ “ 0 the linearized Faraday's law expresses 

— j V . Kjx ttAt;; ci&i-itiatory xiiagiietic iLeid in 

Ic'rms of the nonzero electric field components 


“ «}(*^z^y * ^y^z) 

Substituting (14) into the linearized Ampere's law yields the 
oscillatory current components 


(14) 



(15a) 


(15b) 



We substitute for v^.» n and from (8ab), (9) and (10c), making use 

of (Ha) to eliminate T /T . Equations (15a, b) reduce to a 2x2 matrix (16) 
« e e 

in E and E . C. = B /4nnM is the Alfven speed. The determinant of 
y z A ^ 

(16) must vanish for a nontrivial solution; the cross term cancels much of 
the main diagonal, yielding the collisional electromagnetic dispersion 
relation. This reduces to the fluid limit of the collisionless electromag- 
netic dispersion relation obtained by Coroniti and Kennel (1970) from 
kinetic theory. We will now examine tlie structure of (16) in the large 
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2 2 2 

Electrostatic Limit “* “ 

In the very low 3 <m/M limit appropriate for ionospheric application 
where C <a^ < we let - ■= in (16), and the Alfv^n terms are 

clin.inatcd from the electromagnetic dispersion relation, which reduces to 

i Vj^/b ) + g/Lj^^ 




, [w + i vii ( V " ^ )1 

■ tuj + ivVr 


, 1 + pb] - «uq + I Vj^D 


= 1 + R- 


(17a) 


The second term on the left hand side of (17a) (an electromagnetic correction) 

is neolieible for v »k^C^/(s or o)/ »m/2M. The rest of (17a) corresponds 

y Z S ^ 

to the electrostatic limit which can be obtained directly by assuming quasi- 
neutrality, n. =*n^, and substituting and e!^ = - in (9) and (12a). 

As required by the continuity equation (1) this electrostatic result can also be 
obtained by making the same substitutions in the conductivity matrix J = £ . E 
(16) and setting 7 . J = 0. recalling that = 0. Neglecting the energy 
transfer terms in yand 1, neutral collisions, acoustic terms and gravity, 
this electrostatic dispersion relation (I7a) reduces to that obtained by Tsai 
et al . (1970). 

We have written (I7a) in a form which separates the Interchange 
mode on tlie left liand side from the product of entropy and drift mode terms 
respectively on the right liand side. Note tlmt U.e interchange mode apn.vnr 
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only as FLR and acoustic corrections in the large limit. Also 

note that in the large or small v limit where (X- X, the entropy 

II 

mode (Hudson and Kennel, 1974c) is eliminated from (I7a) which becomes 

(u>+ 

« - lvii{w [1 + pbl - uiQ + ivj.p - (17b) 

Neglecting neutral collisions, acoustic terms, and gravity this quadratic 
isothermal dispersion relation reduces to that obtained by Chu et al. (1969). 

Electromagnetic Flute Mode Limit ^ 

The limit k -* 0 reouires that v -* 0 and (16) becomes 
z ■ II 

I ^ UJ + )((U+lv/b)+g/Lj^| 

{((U) (0)+ (lU- kyVj5)| = 0 (Ig) 

where v X “• 2 (m/M) as v -• 0. 

II ® II 

The first bracket in (18) contains interchange mode terms. The 
second bracket in (18) contains a zero frequency ion acoustic mode (tu), 
a damped entropy mode ( u> + iv ^), and a purely oscillatory drift mode 
( ( 1 ) “ k V ). These three modes are basically electrostatic, and require finite 

' y y 

k for growth. In low 0 <<1 plasmas the Alfven speed >>C^, the 

ion acoustic speed, which is the characteristic electrostatic phase velocity 
alonp the magnetic field. Hence these finite electrostatic modes are 
properly treated in the large k^C_^/cu limit of the electromagnetic disper- 
sion relation {16). 
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2 2 2 

■'J)<;cou plt'<l intcrcliani; e mode: k C./u) <<1 

In t!if pri'ccding section we saw that the interchange mode was 
2 2 2 

ul)laim'fl in tiu* smeill C./u) limit of the coUisional electromagnetic 

Zi A 

2 2 2 

dispe r relation (16). Retaining terms of order k C. /(« and dropping 

2 Z 2 

tnrms of order bk C./ud , the interchange mode dispersion relation is 

Z x\ 


^u) + R|UJq^ (u) + I )+g/Lj^-k^C^ = 0 


(19a) 


The solutions are 
U) 


= l/2| - (RjUJq + I vj_/b) 


+ 4(|<2c 2 - g/L^)] (19b) 

2 2 

’.i’lic k_ C ^ Leriii in {19b) is strongly stabilizing in most applications. 

Solving the electrostatic dispersion relation (17a) for the finite k 

correction to the interchange mode gives a less restrictive and incor- 

2 2 

red rcjsult. The k C. term . rest r lets the interchange instability to 

Z A> 

0 in slab geometry, and to the lowest order flute perturbations of 
an entire flux tube in general geometry, Haerendel (1974) has performed 
a dipole geometry calculation treating magnetic field lines as equipoten- 
tials. Assuming that the plasma density is field-aligned, he has averaged 
all density-dependent quantities in the local dispersion relation (20d) 
over an entire flux tube. Although the plasma density is altitude - 
rather than field - aligned below the F maximum, the step to averaging 
over altitude-aligned plasma density in dipole geometry is not straight- 
forward. We will not pursue the general geometry flute mode further; 
instead we will focus on new results in slab geometry. 
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Mode; = 0 in Slab Geometry 

Setting k =0 and neglecting gravity in (19b), 
and tlie otlier is marginally stable 

u)_ = “ 


one root is damped 


(20a) 


Including gravity, the roots are 


““ “ l{ ‘ (^T®0 ^ 

iAVo- Ivj,/b)2 . 4g/L^} 

ui| becomes unstable and can be written in the form 


(20b) 
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(20c) 


whore 

X - (Rt“d)^ ■ - 4g/Lj^ 


y - 2(Ry»pVji/b) < ,u 

n-OxaO n=»x<0 

At b = 0 the real frequency vanishes and the growtii rate is 
maximum 


u) 




tj! 


I"- 


I n 


■“IV 


Z z 


V t n 
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(20d) 
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In tins limit v ./b - and = 0» bonce the nrowth rate is 

in i u 

independent of Coulomb collisions and FLR effects. This limit corres- 
ponds to that obtained by Haerendel (Balsley et al, 1972; Haerendel 
1971) for slab geometry. We have plotted this maximum growth rate as 
y/ /g / R jL vs from (20d) in Figure 1. We see that the col- 

li sionless growth rate is maximum and that the growth rate decreases 
with increasing neutral collision frequency. 

In the collisionless limit Vjl -• 0, (20c) yields precisely the resott 
obtained by Rosenbluth et al . ( 19 ^ 2 ) from kinetic theory 


where FL.R effects are stabilizing for 


(21a) 




(21b) 

This sets an upper limit on unstable b or lower limit on unstable perpendicular 
wavelengths. 


In the collisional limit of (20c) 




we can expand the radical to obtain one growing root 


V + 1 


(22a) 


(22b) 
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In tlio limit Vj^/b the growth rate is 


Y = I 


_^b_ Vt/b « V f _ + 0. 3bv . . 


(22c) 


wliich agrees with the collisional Rayleigh- Taylor mode growth rate of 
Ilaerendcl (1974; Balsley et al. 1972), here modified to include ion 
viscous damping, which we have seen vanishes for b = 0. 

'i'he full dispersiori relation (20c) is more general tiiari elfu-ir 
i . e result of Haerendel (20d) or Roscr bluth et al . (21a). In adcU-.ion to 
including Coulomb collisions, (20c) demonstrates the effect of neutral 
collisions as well on the collisionless FLR cutoff. In Figure 2 we 
have plotted the complete Rayleigh-Taylor growth rate (20c) at k^= 0 
as a function of b at fixed v.., L and two different values of 

The top curve is essentially collisionless (Vj^ <'2/g/i-»x )J the bottom curve 
is collision-dominated (Vj^ >2/g/L ). The growth rate is maximum at 

b = 0 in both cases, and smaller in the collision-dominated case, as 
indicated in Figure 1. Collision dominated growth extends to shorter 
perpendicular wavelengths than permitted by the collisionless FLR cut- 
off (Rosenbluth et al. 1962) of the top curve. The y term in (20c), which 
vanishes in the coUisionless limit and becomes significant as b -» 1 in the 
collisional case, is responsible for the extension of the instability region 
dt»wn to the ion Larmor radius, where the present fluid analysis breaks 
down. Physically the collisional particle diffusion drift opposes and 
i;anc els the coUisionless FLR drift, Tliis diffusion drift depends on 

v,/b = V -f 0. 3 b V.., hence both neutral and Coulomb collisions at 
' in 11 

finite perpendicular wavelengths. 
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Ion viscosity only affects transport when there exists a spatial 

gradient. For 1 the important spatial gradient is the perpendicular 

wavelength or b. It is clear from {22c) that in a plasma where v v. , 

ii 

such as the nighttime equatorial F region, ion viscosity previously neg- 
lected will dominate ion-neutral collisions in determining the perpendicular 
ion diffusion rate Vj^/b and Rayleigh-Taylor growth rate at short perpen- 
dicular wavelengths. 


Conclusion 

From two-fluid theory we have derived a collisional electromagnetic 
dispersion relation for the interchange mode coupled to the Alfven, acoustic, 
drift and entropy modes. We have demonstrated the fundamental electro- 


magnetic nature of the interchange mode which appears in the opposite 
2 2 

' 2 “A’ ^ ^ X ^ 0 - 1 . c*.^.-Ccri^rnagiieti^^ ciiopersion relation (lo) irom 

the electrostatic modes (k^C^/cv^ >> 1). The main stabilizing affect 
for finite comes from coupling to the electromagnetic Alfven mode 
rather than electrostatic modes, even though interchange mode terms 
appear in the electrostatic dispersion relation. Coupling to the Alfven 
mode restricts the interchange mode to 0 in slab geometry. 

In all limits the interchange mode destabilized by gravity has a 
positive growth rate only when gravity is directed antiparallel to the 
density gradient. For example, this precludes the growth of the Rayleigh- 
Taylor instability on the topside of the F-layer. The growth rate maxi- 
mizes at = 0 and b = 0. Only neutral collisions affect the growth 
rale at b =: 0, and the growth rate decreases with v.^. For finite b, 
ion-ion collisions contribute to the perpendicular ion diffusion rate 

~ ^hi * ^* 3 1) V... 'I’here is competition between tlio collisional 
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,di/ fusion drift and the coil is ionless FL-H drift at short perpendicular 
wavelengths. When < 2 /gi L j_ , the coUisionlcss FLR drift mode 

(Rosciibl iith et al, 19^2) stabilizes short perpendicular wavelengths. 

Wl^en V |^/b >2 /g/F ^ colli sions eliminate the FLR cutoff and extend 
the unstable spectrum down to the ion Larmor radius where the small 
b expansion breaks down. 

Since v,. > > v. at the altitudes where spread F is generally 
11 in L o j 

observed, including ion viscosity is an important extension of the col* 

lisional Rayleigh-Taylor instability theory (Balsley et al. 1972; Haerendel 
1974), Our consideration of finite perpendicular wavelengths has shown 
that it is possible for the Rayleigh-Taylor instability to cover the entire 
range of perpendicular wavelengths above the ion Larmor radius in the 
limit Vj^/b >2 /g/L where collisions eliminate the FLR cutoff. How- 
ever, typically Vj^/b < 2 /g / L ^ at the altitudes where spread F is 
observed (Hudson and Kennel 1974b). Hence, the Rayleigh-Taylor instability 
is primarily coUisionless and Haerendel's limit does not apply. The 
FLR cutoff then restricts the Rayleigh-Taylor instability to perpendicular 
wavelengths the order of a hundred meters or greater, and the drift 

mode which is investigated in the companion paper (Hudson and Kennel 

1974b) has a larger growth rate at shorter perpendicular wavelengths. 
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Table 1 Notation 


1 


Rt = t./t^ 


m 


/M 


V . - 

ei 


(A/10)n^ 

3 72 


3. 5 X 10 T 


(eV) 


parameter perpendicular to 
magnetic field B 

parameter parallel to magnetic 
field B 

perturbed quantities which vary 

i (k . r - lot) 
as e — 

ion electron temperature ratio 

electron- Lon mass ratio 

electron- ion collision frequency, 
single charged ions 


V.. = 'f /2 M ) V . 

11 ei 

V = 3. 5 X lO’^n T (eV) 
en n c 


V. = /(m/2M) V 
in on 


V 

e 


V , + V 

ei en 


A 

a = / 2T /m 

e V e 

A. = J2T./M 
1 V i 



ion-ion collision frequency 

electron- neutral collision 
frequency (Nicolet 1963) 
ion-neutral collision frequency 

total electron collision frequency 

Coulomb Log (Spitzer 1967) 
electron thermal speed 

ion thermal speed 

ion acoustic speed 



0 = eB/Mc 
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ion cyclotron frequency 
electron cyclotron frequency 
ion Larmor radius 


0 = eB/mc 

e 

p. =A./Q 
1 1 


p = a /Q 
e e e 


(U = 4TTne / M 
P 

X = a / V 
e e e 


X.. = A./v.. 
11 1 11 


= (Vn/n) 


-1 


b = k^C^/Q^ 

J’ s 

= -^cT^/eB^^kj^. (V n X B)/n 


V , = b V . +0. 3b V.. 

in ii 

V|i = k^a^/ 2 C V 
" z e re 


C = a, /a 

r 



C 

X 


electron Larmor radius 


ion plasma frequency 

total electron mean free path 

ion- ion mean free path 

perpendicular density gradient 
scale length 

ion finite Larmor radius parameter 


Diamagnetic drift frequency, de- 
fined to be positive for waves 
propagating in the electron drift 
direction 

perpendicular ion diffusion rate 

parallel electron streaming rate 

ratio of perpendicular to parallel 
electrical conductivities (Braginskii 
1966) 

dimensionless thermoelectric 
coefficient 

dimensionless electron thermal 
conductivity coefficient 
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X = (2/3) [c^C^ + (l + cj] 

+ 12 (m/ M) V II 

+ i (l - R^)(m/M) (v^„- 3v^.) 

T= (2/3) (l + gJ^ 

+ z( 1 - Cj,)( 1 - R^ )(m/M) 


heat conciluction and energy 
transfer coefficients 


All other notation standard 
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Figure 1 

1 - 

Plot of maximum Rayleigh-Taylor growth rate in slab geometry 
(20d) which occurs for b = 0 and k =0, normalized as y! / ^ vs. 

v.^/v^ g/ F jL » Note that the growth rate is maximum in the collisionless 
limit and decreases with increasing neutral collision frequency; ion 
viscosity and FLR stabilization do not affect the b = 0 mode. 

We have scaled the x = 'f % axis with altitude for an average 

model neutral atmosphere (Johnson 1965) and the y = y/ / g/Lx axis with 

growth rate for two reasonable choices of Lj_for the bottom of the nighttime 
F layer, = 9. 42 and 25. 6 km. The altitude dependence is affected more 
than the magnitude of the growth rate by this change in L»x» 

Figure 2 

pi rrb-Xaylor gcowth rate (20c) at k =0 as a function of h at 

z 

fixed L 1 = 9. 42 km and v, . =0. 55 sec for two different values of 
^ 11 

v.^ = 0. 014 and 0. 08 sec . The top curve is essentially collisionless 

<7./ g/ the bottom curve is collision-dominated > 2/ g7 L x J * 

The growth rate is maximum at b = 0 in both cases and smaller in the 
collision-dominated limit. Collision dominated growth extends to shorter 
perpendicular wavelengths than permitted by the coUisionless FLR cutoff 
(Rosenbluth et al. 1962) of the top curve, but the theory breaks down as 
b -• 1 (dashed line). 

The two different neutral collision frequencies 014 and 

0. 08 sec correspond to two different altitudes h = 400 and 280 km for 
an average model neutral atmosphere (Johnson 1965). 



Altitude (km) 
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